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Abstract 



We develop a simple theory of the electromagnetic response of a d- wave supercon- 
ductor in the presence of potential scatterers of arbitrary s-wave scattering strength and 
inelastic scattering by antiferromagnetic spin fluctuations. In the clean London limit, the 
conductivity of such a system may be expressed in " Drude" form, in terms of a frequency- 
averaged relaxation time. We compare predictions of the theory with recent data on YBCO 
and BSSCO crystals and on YBCO films. While fits to penetration depth measurements 
are promising, the low temperature behavior of the measured microwave conductivity ap- 
pears to be in disagreement with our results. We discuss implications for d-wave pairing 
scenarios in the cuprate superconductors. 



I. Introduction 



A remarkable series of recent microwave experiments on high quality single crystals 
of YBCO 1 ' 2 '^' 4 '^ has been taken as evidence for d-wave pairing in the high-T c oxide su- 
perconductors, complementing NMR, 6 photoemission, 7 and SQUID phase coherence data 8 
supporting the same conclusion. 9 In particular, there is thus far no alternate explanation 
for the observation of a term linear in temperature in the YBCO penetration depth, 1 
other than an unconventional order parameter with lines of nodes on the Fermi surface. 
Several initial questions regarding discrepancies between this result and previous similar 
measurements, which reported a quadratic variation in temperature, have been plausibly 
addressed by analyses of the effect of disorder, which have suggested that strong scattering 
by defects in the dirtier samples can account for these differences. 10 ' 11 

We have recently attempted to analyze the dissipative part of the electromagnetic 
response, i.e. the microwave conductivity cr, within the same model of d-wave supercon- 
ductivity plus strong elastic scattering, to check the consistency of this appealingly simple 
picture. 12 We found that the conductivity could be represented in a Drude-like form in 
which the normal quasiparticle fluid density and an average over an energy dependent 
quasiparticle lifetime entered. For microwave frequencies small compared to the average 
relaxation rate, the conductivity was found to vary as T 2 at low temperatures approaching 
ne 2 /ttAqui at zero temperature. Here Aq is the gap maximum over the Fermi surface. 
At higher microwave frequencies, the interplay between the microwave frequency and the 
quasiparticle lifetime was found to lead to a nearly linear T dependence over a range of 
temperatures. While some of the qualitative predictions of this model are in agreement 
with experiment, the low-temperature T predictions for the low- frequency microwave 
conductivity differ from the linear-T dependence reported. 

The main purpose of this paper is to explore further the overall consistency of the d- 
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wave pairing plus resonant scattering model predictions for the low-temperature behavior 
of the electromagnetic response of the superconducting state. We will also examine the 
electromagnetic response over a wider temperature regime by phenomenologically including 
the effects of inelastic spin-fluctuation scattering. In the process we intend to provide the 
derivations of results reported in our previous short communication, 12 and address various 
questions raised by it: 

1) To what extent can the microwave conductivity in a d x 2_ y 2-wave superconducting state 
be thought of in direct analogy to transport in a weakly interacting fermion gas with 
a normal quasiparticle fluid density n qp (T) and a relaxation time t(u) characteristic 
of nodal quasiparticles? 

2) Can the temperature dependence of the microwave conductivity be used to extract 
information on the quasiparticle lifetime? 

3) What is the characteristic low-temperature dependence of the quasiparticle lifetime for 
resonant impurity scattering in a d x 2_ y 2 superconductor and how does it affect a? 

4) What happens at higher temperatures when inelastic processes enter? 

5) What happens to &\{T, ft), A(T, ft) and the surface resistance R s (T,ft) at higher mi- 
crowave frequencies? 

6) To what extent can a model with a d x 2_ y 2 gap plus scattering describe the observed 
penetration depth and conductivity of the cuprates? Can the response of a d x 2_ y 2-wave 
state be distinguished from that of a highly anisotropic s-wave state? 

The plan of this work is as follows. In section II, we derive the expressions necessary 
for the analysis of the conductivity and penetration depth of a superconductor in the 
presence of impurities of arbitrary strength within BCS theory. In section III, we examine 
several useful limiting cases of these results analytically. In Sec. IV, we introduce a natural 
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definition of the quasiparticle lifetime which allows the conductivity to be cast in a "Drude- 
like" form with a temperature dependent carrier concentration n qp (T). Then we describe 
results obtained from a model for inelastic scattering by antiferromagnetic spin fluctuations 
and include these in a phenomenological way so as to describe the conductivity over a 
wider temperature regime. In section V, we compare results for the penetration depth, 
conductivity and surface impedance with data on high-quality samples, including both i) 
scaling tests of the d-wave plus resonant scattering theory at low temperatures, and ii) fits 
over the entire temperature range. In section VI we present our conclusions concerning 
the validity of the model and suggestions for future work. 

II. Electromagnetic response: formalism 

We first review the theory of the current response of a superconductor with general 
order parameter to an external electromagnetic field, with collisions due to elastic 
impurity scattering included at the t-matrix level. 13 ' 14 ' 15 We expect such a theory to be 
valid at low temperatures in the superconducting state, if inelastic contributions to the 
scattering rate fall off sufficiently rapidly with decreasing temperature. This is the case in 
the model we discuss most thoroughly, namely a d x i_yi state with an electronic pairing 
mechanism. In such the gap opens, the low-frequency spectral weight of the 

interaction is supressed and the dynamic quasiparticle scattering decreases. The scattering 
rate in the superconducting state contains two factors of reduced temperature T/T c for 
electron-electron scattering, and one for the available density of states in the d-wave state, 
and therefore varies as (T/T c ) 3 at low temperatures. At temperatures of order .3 — AT C 
the dynamic scattering has decreased by one or two orders of magnitude from its normal 
state value, at which point elastic impurity scattering dominates the transport. In this 
low temperature region, the gap is well formed and its frequency dependence occurs on 
scales larger than T c . Thus it is appropriate to model this system within a BCS framework. 
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Furthermore, since the dominant quasiparticle density is associated with the nodal regions, 
we assume that the qualitative features of the temperature dependence of the transport 
will be unaffected by the details of the band structure, and consider a cylindrical Fermi 
surface with density of states Nq, and an order parameter A^ = Ao(T)cos2</> confined 
to within a BCS cutoff of this surface. A more complete theory capable of describing 
the higher temperature regime where inelastic scattering processes become important is 
discussed in section IV. 

If an electromagnetic wave of frequency O is normally incident on a plane supercon- 
ducting surface, the current response may be written 

2 

j(q, O) = -K(q, O) A(q, Q) = - [£„(q, O) - — 1 A(q, O), (1) 

L mc J 

where A is the applied vector potential. The response function is related simply to the 
retarded current- current correlation function, with 

K p (d,n) =<[j,j] i? >(q,0)^ 
- ( ) <kk d£ k T^tr\g(k + ,u n )g(k-,uj n - fl m ) 



- > ^.uf2 m ^il+i0+' 



n 



where k± = k±q/2 and u n = (2n + 1)ttT and Vt m = 2mirT are the usual Matsubara 
frequencies. The approximate equality in the last step above corresponds to the neglect of 
vertex corrections due to impurity scattering and order parameter collective modes. The 
former vanish identically at q = for a singlet gap and s-wave impurity scattering, 16 while 
the latter are irrelevant if the order parameter corresponds to a nondegenerate represen- 
tation of the point group. As usual, in the last step we have performed the analytical 
continuation iO m — > O + z0 + . The single particle matrix propagator g is given as, e.g., in 
Ref. 16 in terms of its components in particle-hole space 

9 (k, u n ) = — T — 3 

W 2 + ^2 + | Afc |2 
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where the r z are the Pauli matrices and Aj. is a unitary order parameter in particle- hole and 
spin space. The renormalized quantities are given by ui n = u n — Eo(u; n ), £fc = £fc + ^(^n), 
and Afc = A^ + Ei(u; n ), where the self-energy due to s-wave impurity scattering has been 
expanded E = E^r*. The renormalization of the single-particle energies ^ measured rel- 
ative to the Fermi level is required for consistency even in the s-wave case, although it 
is frequently neglected because in the Born approximation for impurity scattering such 
renormalizations amount to a chemical potential shift. For a particle-hole symmetric sys- 
tem, these corrections can be important for arbitrary scattering strengths, but are small 
in either the weak or strong scattering limit. 16 ' 17 We therefore neglect them in what fol- 
lows, and postpone discussion of the particle-hole asymmetric case, where these effects can 
become large, to a later work. 

A further simplification arises for odd-parity states and certain d-wave states of current 
interest, where a reflection or other symmetry of the order parameter leads to the vanishing 
of the off-diagonal self-energy Ei. In this case, the gap is unrenormalized (Aj. = A&), 
leading to a breakdown of Anderson's theorem and the insensitivity of the angular (e.g., 
nodal) structure of the gap to pairbreaking effects. 

Rather than solve the self-consistent problem in full generality, in most of what follows, 
we focus on two cases of special interest: i) s-wave pairing with weak scattering, for 
purposes of comparison; and ii) d-wave pairing without Aj~ renormalization for weak or 
resonant s-wave scattering. In case i), the self-energies Eq = T^Gq and Ei = Tj^G\ 
are the familiar integrated Green's functions from Abrikosov-Gor'kov theory, where r^y 
is the scattering rate at T c attributable to impurities alone, and we have defined G a = 
(i/27rNQ)T,i i Tr[T a g]. The Green's function (3) and the self-energies must be calculated 
together with the gap equation, A(k) = T J2 n Sfc' ^ / fcA:'T r ( r l/2)fi , (^ / , w n ), where V^i is 
the pair potential. In Sees. II-III, all calculations are done self-consistently within weak- 
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coupling BCS theory, which yields Aq/T c = 2.14 for a pure d x 2_ y 2 state. When comparing 
with experimental data in Sees. IV-V, we adopt larger values of Aq/T c of 3 or 4 to simulate 
strong-coupling corrections. 

We now continue the derivation of the response on a level sufficiently general to sub- 
sume both cases i) and ii) above. If we neglect ^ renormalizations, the self-energies are 
given in a t-matrix approximation by 

_ TG _ -Tgi 

E ° " c^ + ^-Go 2 ' 1 " c2 + Gl 2 -G 2 ' (4) 

where V = ninj (itNo) is a scattering rate depending only on the concentration of defects 
the electron density n, and the density of states at the Fermi level, Nq, while the strength 
of an individual scattering event is characterized by the cotangent of the scattering phase 
shift, c. The Born limit corresponds to c > 1, so that T/c 2 ~ F N , while the unitarity 
limit corresponds to c = 0. To evaluate Eq. (2), we first perform the frequency sums, then 
perform the energy integrations as in Ref. 15, yielding in the general case 



Reif(q, Q) 



1 ne z 



2 mc 



— k : k / du{ [tanh- tanh/3 

2tt J U 2 2 

+ tanh— + tanh-^ 

L 2 2 



Re/-| — (a;, u — 0,) + 
Re/++(a;, u — 



(5) 



Imif(q, O) = - 



1 ne z 



2 mc / 2-7T 



- k : k du 



tanh- tanh-^ x 

2 2 J 

xIm{f-)-+(a;, cj — O) — 7-| — (a;, a; — 



(6) 



In calculating the surface impedance of the cuprate superconductors, it is important 
to take into account the anisotropy of these layered materials. 18 Here we are interested 
in the response associated with currents which flow in the ab layers. The wavevector in 
the ab plane is determined by the long wavelength of the microwaves and hence can be 
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set to zero. Furthermore, the short quasiparticle mean free path in the c-direction means 
that the surface impedance is determined by the conductivity of a CuC>2 layer. Thus the 
surface impedance in this case is given by 

z4ttO \ 1/2 



\c z {ai{n,T) -ia 2 {yt,T)J 

Here a\ — ia 2 is the complex frequency- and temperature-dependent q = layer conductiv- 
ity. It is customary to write the imaginary part of the conductivity in terms of a frequency- 
and temperature-dependent inductive skin depth A(f2, T), 

c 2 

° 2 = 4^(0, TY (8) 

At temperatures a few degrees below T c , a 2 3> <j\, so that the surface resistance R s is 
given by 

R s = ReZ(n,T)^ W(n -P^ n - T \ (9) 

c 4 

and the surface reactance X s is 

X, = Im Z (fl,r)^ 4 ^ A f T > . (10) 

Thus microwave surface impedance measurements provide information on the inductive 
skin depth A(fi, T) and the real part of the conductivity ai(0, T). In the previous section, 
we have dropped the subscript 1 and denoted the real part of the conductivity simply by 
cr(fi, T), and in the limit O — > 0, A(0, T) is just the London penetration depth. 
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At Q = 0, the energy-integrated bubbles I++ and 7-| are given by 



j ++(cj J) = J_ _ ^(^+ + ^+)+A / fc+ (A fc+ -A / fc+ ) 
£o+ (Co • • iii . )£o • Co . 



and 



(11) 
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/+-(u/,u/) = — + *-(*++*-)+ A U A fc+- A U (12) 

£0+ (£0+ -£o-)£o+£o- 

Here u a = a)(a; + za0 + ), A^q, = A^.(a; + ia0 + ), and £oa = s g n 0J \j^ J 'a ~ with a = ±1. 

We first consider the dissipative part of the response, reflected in the q = conductivity 
a(O) = -(c/n)ImK(q = 0,O). Combining Eqs. (6,11-12) yields 



ne 2 /" 00 1 1 

^j(fi) = y ^{tanh[-^] - tanh[-/?(u; - CI)]} x S^oj, O), (13) 



where 



Sij(u, O) = im 



2tt 1 J 



^((^ + + ^)+A / fc+ (A fc+ -A / fc+ ) , i i 



(e 2 +-ef+) W 

u'_(Q+ + u'_)+A' k _(A k+ -A' k _) f i i 



+ 



«» 2 + 



£()+ fo-' 



(14) 



and primed quantities are evaluated at u> — O. For <i-wave pairing there is no gap renor- 
malization, so that Aj~ a = A^ and the kernel Sij reduces to 



Sjj(u),Q) = Im J —kikj x 



1 1 



+ — 



to. 



1 1 

+ 



• (15) 



We also require an appropriate expression for the London limit Meissner kernel 
Re K(0, 0) to evaluate the penetration depth. Taking O — > in Equation (5), we obtain 19 ' 20 



ne 



mc 



dw tanh 



2n 3 \£3 + J 



(16) 



In the special case of isotropic s-wave pairing and Born scattering this reduces to the well 
known result 21 ' 22 

Re K '(0,0) = -— / duo tanh^ Re( — ^ = ), (17) 

K ' 2mcJ 2 l( v 2_ A 2 )[Vu 2 - A 2 + zI>]J V ' 

with d = (D+A/A. 



III. Limiting cases 

We are primarily interested in the low-temperature, low-frequency conductivity re- 
quired to discuss experiments in the microwave regime. Since the microwave energy is 
generally lower than the temperatures of interest, it is useful to replace (tanh/3u;/2 — 
tanh/3(co> — fi)/2)/(2fi) by its small fi/T limit —df/duo, providing an exponential cut- 
off above the temperature T in the integral (12). At low temperatures T <C Aq, the 
temperature dependence of the conductivity depends strongly on the lifetime of the low- 
energy quasiparticle states, determined by the self-consistent solution to Co = uo — Eg and 
Afc = Afc — Ei, where Eq and Ei are given by Eq. (4). 

In an ordinary superconductor with weak scattering, only the exponentially small num- 
ber of quasiparticles above the gap edge contribute to absorption. Resonant scattering, 
such as occurs in the case of a Kondo impurity in a superconductor, is known to give rise 
to bound states near the Fermi level, reflected in a finite density of states at uo = and 

no 

leading to absorption below the gap edge. A similar phenomenon occurs in unconven- 
tional superconductors, with the difference that, whereas in the s-wave (Kondo) case the 
bound state "impurity band" is isolated from the quasiparticle density of states above the 
gap edge, in unconventional states with nodes the "bound state" lies in a continuum, and 
the lifetimes of all states are finite. 24 ' 2 ^. Nevertheless the energy range between zero and 
the gap edge Aq may be partitioned crudely into two regimes, separated by a crossover 
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energy or temperature T* dependent on the impurity concentration and phase shift. Be- 
low uj ~ T*, the scattering rate — 2ImT,Q(uj) is large compared to u, and the effects of 
self-consistency are important. The physics of this regime is similar to gapless supercon- 
ductivity as described by the well-known Abrikosov- Gor'kov 26 theory of pairbreaking by 
magnetic impurities in ordinary superconductors. The low-temperature thermodynamic 
and transport properties are given by expressions similar to analogous normal state expres- 
sions, with the usual Fermi surface density of states Nq replaced by a residual density of 
quasiparticle states no = N(uj — > 0) in the superconductor. Above T*, self-consistency can 
be neglected, and transport coefficients are typically given by power laws in temperature 
reflecting the nodal structure of the order parameter. 27 We note that this "pure" regime 
will correspond to the entire temperature range if the impurity concentration is so small 
that T* -> 0. 

In this paper we focus primarily on the case of resonant scattering in an attempt to 
describe the physics of Zn doping in the cuprate superconductors. While Zn impurities 
are believed to have no, or very small, magnetic moments 28 , they nevertheless appear to 
act as strong pairbreakers. 28,6 A possible explanation for this strong scattering could be 
associated with the fact that an inert site changes the local spin correlations of its nearest 
and next nearest neighbors 2 ^ These changes can lead to strong scattering^ and even to 
bound state formation 31 for the holes of the doped system. With this in mind, here we 
assume that a Zn impurity may be approximated by an isotropic potential scatterer with 
a large phase shift close to tv/2. 

The essential physics of gapless transport in unconventional superconductors was dis- 
cussed in the context of heavy fermion superconductivity by Hirschfeld et al. 24 and Schmitt- 
Rink et al. 25 Although both works presented calculations for model p-wave states, most 
conclusions reached regarding p-w&ve states with lines of nodes continue to hold for the 
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d-wave states in quasi-two-dimensional materials of interest here. For example, the normal- 
ized density of states N(oj) = — ImGo(cu) is linear in energy for the pure system, and varies 
as no + aT 2 for T <C T* for an infinitesimal concentration of impurities. Neresesyan et 
al. 32 have recently called into question the existence of the residual density of states ng in 
a strictly 2D system. We believe nevertheless that both the underlying three-dimensional 
character of the layered cuprates, as well as the extremely low temperature at which the 
difference between the logarithmic term and the slow power law behavior found in Ref. 32 
becomes significant, make such considerations irrelevant for our purposes. 

All quantities of interest in the gapless regime may be obtained by expanding Co (and 
A if necessary) for u < T*, with the result Co ~ 1(7 + 6a; 2 ) + au, where 7, a, and b are 
constants. T* itself may be shown to be of order 7. In the case of a d x i_ y i state over a 
cylindrical Fermi surface, 7 satisfies the self-consistency relation 7 = rno/(c 2 + no 2 ), where 
no = 2//tK(zAq/7), with K the complete elliptic integral of the first kind. For small im- 
purity concentrations such that V <C Aq, one finds no — (27/71 Ao) ln(4Ao/7). In the Born 
limit, c ^> 1, 7 ~ Tj^uq, an d both 7 and no therefore vary as ~ Ao exp(— Aq/Tj^). In the 
resonant scattering case of primary interest, on the other hand, 7 = T/uq and for small con- 
centrations the residual scattering rate is determined by (7/A0) 2 = (7iT)/[2Aorn(4Ao/7)]. 
The constants a and b are found to be \ and —1/(87), respectively. Thus for strong scat- 
tering both 7 and the residual density of states uq vary as (TAq) 1 / 2 up to a logarithmic 
correction. This is important because it means that low-energy states may be strongly 
modified, even though the impurity scattering rate, which varies as V near T c , is insuf- 
ficient to suppress T c significantly. In the usual Born limit, on the other hand, gapless 
effects become important only when ~ Aq, implying a large T c supression. As the 
normal state inelastic scattering rate, of order T c in temperature units, is much larger 
than the impurity scattering rate in clean samples, we expect that impurities are in any 
case relatively ineffective in suppressing T c until the elastic scattering rate at the transition 
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becomes a significant fraction of the inelastic one (see Sections IV and V). 



These estimates enable an immediate evaluation of Eqs. (13) and (15) in the gapless 
regime, 



where ago = ne /(m7rAo(0)) for a d x 2_ y 2 state. The first term in Eq. (18) is a remarkable 
result first pointed out by P. A. Lee, 33 namely that the residual conductivity a(Q — > 0, T — > 
0) of an anisotropic superconductor with line nodes on the Fermi surface is nonzero and 
independent of impurity concentration to leading order. It arises technically from the 
first term on the right hand side of Eq. (14), and is present in principle regardless of the 
scattering strength. Physically this reflects a cancellation between the impurity-induced 
density of states and the impurity quasiparticle scattering lifetime. The linear variation 
w/Aq of the d-wave density of states is cut off when u drops below the impurity scattering 
rate r _1 . Therefore, at low energies there is a finite impurity-induced density of states 
which varies as (Aqt) -1 . At low temperatures such that T < r _1 , the effective relaxation 
rate which determines the conductivity is proportional to the density of states (Aqt) -1 
multiplied by r, giving A^ 1 independent of the scattering strength. Very recently it was 
pointed out that a generalization of the present theory to include a finite scattering range 
results, in the limit of sufficiently large range or disorder, in a residual conductivity which 
scales with the scattering time (2r) -1 . 34 The predicted residual conductivity in this regime 
is however too small to apply to the experiments considered here. 

In Figures 1 and 2 we illustrate the effect of varying the phase shift and impurity 
concentration on the T-dependence of the conductivity with a full self-consistent numer- 
ical evaluation of Eqs. (13) and (15) for a d x 2_ y 2 state. The intrinsic gapless behavior 
represented by Eq. (18) is clearly visible in the resonant limit, c ~ 0, but in the Born 
limit, c ^> 1, the same limiting behavior is effectively unobservable for small concentra- 




(18) 
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tions at O = 0. Instead, the conductivity tends to a value gq = ne 2 /2mY^ except at 
exponentially small temperatures, where it again approaches a 00 , due to the narrow width 
7 ~ Aq exp — Ao/r^Y of the gapless range in this limit. 

For T > T* ~ 7, we take uj — uj ~ £o( w ) rather than Eo(u)), and keep only the 
leading singular terms in Eqs. (13) and (15) as V — > 0, arriving at the remarkably simple 
expression, 

"~< n > ^ (£) SS<^) N ^ HtT^S))- < 19 > 

where t~^(uj) = — 21m Eq(u;), for any choice of phase shift. Note that N(uj) is the density 
of states for a pure superconductor normalized to N(0) and varies as |u;/Ao| for a d x 2_ y 2 
state at low energies. Eq. (19) is exactly the result expected for the conductivity of 
noninteracting fermions with density of states N(uj) and 1-body relaxation time t(u), and 
is reminiscent of the Drude-like expression used by Bonn et al. to analyze their data. 
However, as pointed out in Ref. 12, the ^-dependence of the superconducting density of 
states tends to induce a strong energy dependence in t(uj) in either the strong or weak 
scattering limits. For a d x 2_ y 2 state we find 

-1, s „ f (7r 2 rA )/[2cln 2 (4A /c)] c ~ 

1 J "l (4T N u/irA )\n(4A /u) c>l 1 J 

leading to the pure limit conductivity result for O <C TAq/T, T <C T c , 

f 2 (T )2ln 2 4Ao C ^ Q 

g xx (Q = o,t)~{ 5 <W r (21) 



CTQ C > 1 

In the opposite limit O > TAq/T, T < T c we find 

a ra (fi,T)^ \>J drT2 * • (22) 

l (^)w ln ^ c»l 

It is instructive to compare the form of the previous results with the more familiar 
form of those expected for an s-wave superconductor with weak potential scattering. We 



begin with Eqs. (13) and (14), and proceed as before in the pure regime, neglecting self- 
consistency in Eq and Ei. We find 

ff -^ K (^) 2 r^(^) w(w)to (n^M)' (s - wave ' Bom) (23) 

where now however the quasiparticle relaxation time in the s-wave superconducting state is 
given by (2r) _1 = — ImEo(a>) — (A/u;)ImEi(u;), and N(uj) = uj/Vuj' 2 — A 2 . This relaxation 
rate has a similar form to that found, e.g., , by Kaplan et al. 35 for the electron-phonon 
quasiparticle relaxation in ordinary superconductors. In the limit O — > 0, T — > 0, we find 

which is similar in form to the well-known Mattis and Bardeen result. 36 

The hydrodynamic limit results Eqs. (21) predict a T 2 behavior 12 for resonant scat- 
tering or a constant 37 behavior for weak scattering for the low-T conductivity of a d-wave 
superconductor under the assumptions set down above. Neither of these is consistent with 
the linear-T variation reported in experiment, which would correspond to the assumption 
of a constant relaxation time r. Thus the low-temperature experimental results appear to 
be inconsistent with the simplest <i-wave model. 12 However, different physical relaxation 
mechanisms than those considered here could change the low-temperature behavior. 

The crossover regime between the hydrodynamic (Eq. (21)) and collisionless (Eq. (22)) 
limits is an interesting one which we investigate further here. In Fig. 2, we illustrate this 
crossover in the Born limit for a d x i_ y i gap, demonstrating that the result a xx — > ctq holds 
only in the hydrodynamic regime fi T^y. This is a point of some importance, since 
experiments on Zn-doped samples appear to indicate a residual conductivity a(T — > 0) 
which scales inversely with impurity concentration, reminiscent of the zero- frequency Born 
result Eq. (21). On the other hand, Figure (2) shows that this behavior disappears at 
microwave frequencies comparable to those used in the experiments. It therefore appears 
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unlikely to us that an explanation in terms of weak scattering can be compatible with the 
observations reported in Ref. 4 and Ref. 5. 



In Fig. 3, we plot the low-temperature conductivity for the case of resonant scattering 
to display the same crossover. It is interesting to note that a quasilinear behavior is in 
fact obtained over an intermediate range of temperatures when the frequency becomes 
comparable to the scattering rate, but this behavior does not appear to hold very far from 



To close the discussion of the low-energy behavior of the conductivity, we give analytical 
results for the frequency-dependent conductivity at zero temperature. 38 In this case the 
factor (tanh/3u;/2 — tanh/3(cj — 0)/2) appearing in Eq. (13) reduces to a window function 
limiting the range of integration from to O. The result may be expanded for small values 
of the integration variable, yielding in the resonant limit 



In Fig. 4, we plot the frequency dependence of the T = conductivity in the impurity- 
dominated regime. 

A full analysis of surface impedance measurements requires, in addition to the conduc- 
tivity a, a knowledge of the inductive skin depth A(fi, T), which reduces in the limit O — > 
to the usual London penetration depth A(T). The = penetration depth in a d x 2_ y 2 
state in the presence of resonant impurity scattering has been calculated by several authors. 
In the gapless regime T < T*, the linear-T behavior characteristic of a <i-wave system is 
destroyed, and one finds the result A ~ Aq + 7rAoT 2 /(67Ao), where Aq = \J mc 1 / Aixne 2 ' 
is the pure London depth, and the renormalized zero-T penetration depth is given by 11 



(A - A )/A ~ (7/(ttA )) ln(4A /7) ^ T/(2 7 ). At higher temperatures T* < T < T c , 
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n ~ r. 




(25) 



the penetration depth crosses over to the pure result, A(T) ~ Ao[l + ln2(T/Ao)]. For com- 
pleteness, we show in Figure 5 the increase of the zero-temperature London penetration 
depth for large values of the scattering parameters in the Born and unitary limits. These 
results are in agreement with those of Kim et al. 39 

The presence of low-energy quasiparticles can induce a strong frequency dependence to 
the low-temperature inductive skin depth A(T, O), which can in some cases mimic shifts in 
low-temperature power laws. Some of these effects were explored in the context of heavy 
fermion superconductivity. 40 Here we observe that the skin depth temperature dependence 
can be suppressed if the microwave frequency is large enough such that Or > 1. In this case, 
it is necessary to use the penetration depth measured at O rather than the limiting low 
frequency penetration depth, to extract the conductivity from surface resistance data. A 
simple expression for the frequency-dependent penetration depth A(T, O) may be obtained 
in the pure regime, T > T*, by neglecting self-consistency in the imaginary part of the 
conductivity as well, 

df\ r (Or) 2 



u) 1 ~a 

OUJ 



(26) 



.l + (Or) 2 J. 

In the collisionless limit Or > 1, the response of the system is perfectly diamagnetic in 
this approximation, A(T, O) — > Ao- In Fig. 6, we explicitly illustrate the effect of increasing 
the microwave frequency on the skin depth of a clean d x 2_ y 2 superconductor. 

IV. Spin fluctuation model for quasiparticle relaxation 

As discussed in Sec. Ill, in the "pure" limit where T* <^ T <^ T c , we find a "Drude"-like 
form (19) for the conductivity of a d-wave superconductor, with r _1 (u;) = — 2ImEo(o;) 
and N(u) the superconducting density of states. In this limit the penetration depth for a 
d x 2_yi state is given by 
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Then using (A(0)/A(T)) 2 = 1 — n qp (T)/n to define a normal quasiparticle fluid density, a 
may be written as 

^(0) = ^g^lm( 1 V (28) 
m \S2 — i/t(u>) I 

where the average (...) is defined by 

(A(u)) = — V ■ (29) 

In the limit where Qt(u) <C 1, Eq. (28) reduces to cr^ = n qp {T)e 2 {r) jm. 

For a d x 2_ y 2 gap, n qp (T) varies linearly with temperature at low temperatures. Thus if 
the average lifetime (r) were constant, a X x would vary linearly with T at low temperatures. 
However, the impurity scattering lifetime is frequency-dependent due to the frequency 
dependence of the single-particle density of states. In Fig. 7 we show plots of t~ 1 (u) 
versus u for the case of a d x 2_ y 2 gap and various values of the scattering phase shift. In 
the unitarity limit we have 

I ( 2T* uj <T* 

-r-\ - \ n 2 TA (30) 

r ( UJ ) I s uj>T* 

2u ln 2 (4A /u;) 

Thus in the "gapless" regime, to < T*, the impurity scattering rate saturates at 2T* and 
in the "pure" regime, u > T* , r varies linearly with u> to within logarithmic factors. In 
this limit, as discussed in Sec. II, the conductivity rises with increasing temperature as T 2 
times logarithmic corrections. This type of behavior is characteristic of a d x 2_ y 2 gap and 
resonant impurity scattering. One power of T comes from n qp (T) and the other from (r); 
both ultimately reflect the linear uj variation of the single-particle energy density of states. 

At higher temperatures, inelastic scattering and recombination processes determine the 
quasiparticle lifetime. In models in which the d x 2_ y 2 pairing arises from the exchange of 
antiferromagnetic spin-fluctuations, 41 it is natural to expect that antiferromagnetic spin 
fluctuations rather than phonons provide the dominant inelastic relaxation mechanism. 
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Calculations of the quasiparticle lifetime 42 have been carried out for a two-dimensional 
Hubbard model in which the spin-fluctuation interaction is taken into account by intro- 
ducing an effective interaction 

3 U 



Here U is a renormalized coupling, and 



nq,*) = - jjVcsf — ( 31 ) 



^BCS (n . A 1 f 1 [i , e P+g e P + A p+q A p ] f( E p+q) ~ f( E p) 



P 

+ 



1 

+ 4 



p+q-^p J U ~ (Ep+q ~ E p) + i0 + 

1 L gp+gCp + A p+g Ap i 1 - f(Ep +q ) - f(E p ) ^2) 

_ e p+q cp + A p+g Ap i f(E p+q ) + / (-Ep) ~ 1 ] 
~ J oo - (E p+q + Ep) + i0+ J 



is the BCS susceptibility with .Ep = ye p 2 + A p 2 , where e p = — 2t(cosp x + cosp y ) — p 
With the interaction given by Eq. (31), the lifetime of a quasiparticle of energy u and 
momentum p in a superconductor at temperature T is given to leading order by 

r in 1 (p^) = -^Yl 



N- 

p' 



/ dv ImV(p-p' 1 v)S(oj-v-E pl )(l + P p )( n (v) + l)[l- f(u - v)] + 

l JO V LO{LO — V)J 

+ f dv Im V(p - p', v)b{v -uj- Ep,) (l - ApA P' ) ( n {y) + l)f(v - u)+ 
Juj+\A' p \ v uj[v-u)J 

roc / A A / \ i 

+ / dv Im V(p - p', u)6(u + v-E pl )(l+ P p )n(v)[l - f(u + u)] 
JO v uj(uj + v)/ J 

(33) 

Here n(v) and f(w) are the usual Bose and Fermi factors, and a quasiparticle renormal- 
ization factor has been absorbed into V. The second term of Eq. (32) corresponds to a 
process in which two quasiparticles recombine to form a pair with excess energy emitted 
as a spin fluctuation. The first and third terms describe scattering processes associated 
with the emission or absorption of spin fluctuations, respectively. 

Quinlan et al. 42 numerically evaluated Eq. (31) to obtain the quasiparticle lifetime 
using parameters for U, t, and the band filling which had previously provided a basis for 
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fitting the nuclear relaxation rate of YBCO 43 and gave a normal state quasiparticle lifetime 
t~ 1 (T c ) of order T c . The temperature dependence of the inelastic quasiparticle lifetime for 
a d x i_ y i gap with 2Aq/T c = 6 to 8 was found to be in reasonable agreement with the higher 
temperature transport lifetime determined by Bonn et al. . At reduced temperatures below 
T/T c of order 0.8, the d x 2_ y 2 gap is well-established and the occupied quasiparticle states 
are near the nodes. Setting p to its nodal value and ui = T, Quinlan et al. found that the 
temperature dependence of the numerical calculations of the quasiparticle lifetime varied 
as T 3 , reflecting the available phase space. 

Figure 8 incorporates results for (r) obtained by setting the scattering rate equal to 
the sum of the impurity and inelastic rates. This procedure neglects the real parts of the 
self-energy as well as vertex corrections arising from the dynamic processes. Nevertheless, 
it shows the qualitative behavior of (r) versus T/T c . Combining a simple parameterized 
fit of the numerical results of Ref. 42 for r" 1 (T) with the unitary elastic scattering rate, 
corresponding results for a(T) versus T/T c are shown in Fig. 9. Here the peak in cr(T) 
arises from the rapid drop in the dynamic quasiparticle scattering rate as the gap opens 
below T c and spectral weight is removed from the spin-fluctuations. 44 The low-temperature 
T 2 dependence implies that at these energies, the quasiparticle scattering rate is increasing 
as the temperature is lowered due to the linear decrease in the single-particle density of 
states and the fact that r is proportional to this density of states in the unitary scattering 
limit. 12 As the microwave frequency O is increased, the temperature T p , at which the peak 
in cr(fi, T) occurs, increases. At the same time the peak value decreases. Adding the 
numerical results for the inelastic scattering rate t^{T) to the unitary elastic scattering 
rate and evaluating Eq. (25) for various microwave frequencies, we find that T p /T c and 
cr(fi, Tp)/<r(0, T c ) vary with O as shown in Fig. 10. 
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V. Analysis 



Quantitative comparison of the simple theory presented here with existing data is 
useful but dangerous. We remind the reader that many features of the model are certainly 
oversimplified, including but not limited to the neglect of the real Fermi surface anisotropy, 
higher-order impurity scattering channels, and strong coupling corrections. However, we 
do not expect inclusion of these aspects of the physics to qualitatively alter the nature of 
the temperature power laws in the response functions at low temperatures in the gapless 
and pure regimes. At higher temperatures T < T c , it is natural to expect that real-metals 
effects will produce nonuniversal behavior in the superconducting state even if the normal 
state is a strongly renormalized Fermi liquid. With these remarks in mind, we proceed as 
follows. We first attempt to fix the impurity scattering parameters within the resonant 
scattering model by comparison to the penetration depth data of Bonn et al. 5 on Zn-doped 
samples of YBCO. It turns out the fit obtained is relatively good in this case, although 
the scattering rates in the case of the Zn-doped samples are not fixed with high accuracy 
because of uncertainties in the zero-T penetration depth. As discussed below, a different 
kind of scaling analysis can be performed on the thin film data of Lee et al. 45 

As one knows from the heavy fermion superconductivity problem, claims to determine 
the gap symmetry by fitting a theoretical prediction to a single experiment on a single 
sample should be treated with caution. It is extremely important to correlate results on 
different kinds of measurements on different samples. The results of the British Columbia 
group afford an excellent opportunity to do this kind of cross-checking. We therefore 
adopt for the moment the "best" results for the scattering parameters in the pure and Zn- 
doped samples from the penetration depth analysis, and use them to compare calculated 
conductivities and surface resistances with the Bonn et al. data. The behavior of the 
temperature-dependent conductivity is much richer than that of the London penetration 
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depth, so it will be important for the consistency of the theory to see which aspects can 
be reproduced by the d-wave plus resonant scattering (plus inelastic scattering) model. 

In Fig. 11, we show one possible fit to the UBC penetration depth data. 5 The curves 
represent the theoretical penetration depth A(T) normalized to the pure London depth 
Aq for different values of the resonant scattering parameters V as given. The value 
Aq/T c = 3 is chosen from the fit of the asymptotic pure d x 2_ y 2 penetration depth 
AA(T) ~ Aohi2(T/Ao) to the intermediate linear-T regime in the pure data (symbols). 
The value T/T c = 8 x 1CP 4 is then chosen by fitting the curvature of the T 2 contribution at 
the lowest temperatures. As the absolute scale of the experimental A(T = 0) is uncertain, 
we have chosen to add constant offsets to the various data sets to try to achieve reasonable 
fits. Figure 11 shows that it is possible to find a consistent choice of such offsets, since 
the scattering rates used for the two Zn-doped data sets, T/T c = 0.018 and 0.009 are in 
the ratio 2:1 as are the nominal Zn concentrations 0.31% and 0.15%. However, a roughly 
equally good fit may be obtained using scattering rates of, e.g., T/T c = 0.03 and 0.006, 
which would then not be consistent with the theoretically predicted scaling of V with the 
impurity concentration n^. Clearly there is a relatively large range of acceptable scatter- 
ing rates corresponding to the two Zn-doped curves, possibly a factor of two or more. A 
determination of the zero-temperature limiting penetration depths of pure and Zn- doped 
samples from, e.g., fiSR experiments, is needed to fix these values more precisely or rule 
out such a fit. 

A procedure for fixing the zero-temperature penetration depth relative to the single- 
crystal data without new experiments has been suggested by Lee et al. They assume that 
the data for their YBCO films follow a universal curve given by the form of the single 
crystal penetration depth in the intermediate temperature regime, as suggested by the 
resonant scattering analysis. Using data on several films, they show that such a scaling 
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is indeed possible, and assign zero-temperature penetration depth values to several films 
on this basis. This allows an internal consistency check of the resonant scattering hy- 
pothesis, wherein one may check to see that the measured coefficients of the T 2 term in 
the penetration depth, equal to C2 = 7rAo/(67Ao) for a d x 2_ y 2 state and resonant scat- 
tering, scale appropriately with the zero-temperature penetration depth renormalization, 
(Aq — Aq)/Aq ^ (7/(77 Aq)) ln(4Ao/7) ^ T/{2 r )). Since a given film in the resonant scatter- 
ing limit is characterized simply by its impurity concentration through the parameter 7, 
using the above expressions it is possible to check scaling without knowledge of the actual 
defect concentration. For example, in Fig. 12 we plot (Ao — Ao) vs. l/c2 for two "different" 
films measured in Ref. 45 actually the same film before and after annealing (films A and 
A' of Ref. 45 ). Each cluster of points in Figure 12 represents a single film, the individual 
points corresponding to differing assumptions regarding other constants, such as the abso- 
lute value of the pure penetration depth, which enter such an analysis. It is seen that the 
agreement with the theoretical scaling is remarkably good, and that this agreement is not 
particularly sensitive to varying assumptions on the subsidiary constants. 

Next we explore whether an equally good fit is possible for the resistive part of the 
conductivity which was also measured in Ref. 5. As we have seen, even in the "pure" limit 
T > T* the conductivity depends on the quasiparticle lifetime. At low temperatures, elastic 
scattering from impurities determines this lifetime. At higher temperatures, however, 
inelastic scattering processes become important and we use a simple parameterized fit 
to the numerical results for the inelastic scattering rate r _1 (T) obtained by Quinlan et 
al. 42 As previously discussed, the parameters of the spin-fluctuation interaction used in 
this work were used in fitting the NMR data and the overall strength was adjusted to 
give t~}(T c ) of order T c . The total scattering rate is taken as the sum of the elastic 
and inelastic rates. Using the usual expression, for the surface resistance R s in terms 
of the real part of the conductivity a and the penetration depth, R s = (87r 2 2 A 3 a)/c 4 , 
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Bonn et al. extracted the conductivity for the same samples whose penetration depth is 
plotted in Fig. 11. In Figs. 13 and 14, we show the conductivity plotted for these samples 
calculated using the elastic scattering parameters taken from Fig. 10 and the inelastic 
scattering results from Fig. 8. Although the size, position, and scaling with frequency of 
the prominent maximum in the conductivity are reproduced qualitatively, it is clear that 
the low-temperature behavior of the data does not correspond to the predictions of the 
model. In section II, we pointed out that, while a a ~ T behavior can be obtained in the 
pure regime if Or ~ 1, it is not generic to the theory; by contrast, the data for at least 
the "pure" sample and 0.15% Zn appear to follow a low-temperature linear-T law for all 
the samples shown. A similar behavior is observed in YBCO thin films and BSSCO single 
crystals. 46 

The further difficulty apparent from the data shown in Figs. 13 and 14 is the rather 
large residual value of the conductivity as T — > exhibited by all data sets. While the 
(i-wave theory predicts a residual absorption, the limiting a 00 — ne 2 /m7rAo of the theory 
is an order of magnitude or so lower than that extracted by the British Columbia group. 4 ' 5 
While qualitatively different physical scattering mechanisms than those considered here, or 
a completely different picture for superconductivity in the cuprates might be responsible 
for the deviations from theory apparent in the data, we prefer to reserve judgement until 
further data is available. Very recent results from the British Columbia group indicate that 
twin boundaries may be responsible for the residual conductivities observed, and possibly 
also account for part of the temperature dependence observed at low temperatures. In 
Fig. 15 we show data for a twin-free, high-purity YBCO crystal 5 compared to the same 
theoretical prediction used for the low-frequency conductivity displayed in Fig. 13. It is 
evident that the residual conductivity in the untwinned has been dramatically reduced, 
and the low-temperature fit to the <i-wave theory correspondingly improved. Clearly high- 
quality Zn-doped samples of this type are also desirable. 
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For completeness we also calculate and display the surface resistance R S (T) for var- 
ious values of the scattering parameters in Figure 16. Here again, we see that the low- 
temperature behavior of the theory is in disagreement with the data. 2 ' 3 ' 4 ' 5 This reflects 
the much lower residual conductivity predicted for our model, as well as the T 2 power 
law dependence. In addition, in order to reproduce the dramatic decrease in R s which 
is observed below T c , we need a large Aq/T c = 4 ratio. It is also important in making 
this comparison to recall that the drop in R s just below T c reflects less the collapse of 
the inelastic scattering rate which enters the conductivity a than the divergence of the 
penetration depth depth near T c (recall R s ~ A 3 ). The data suggests that the magnitude 
of the gap opens more rapidly than usual. This type of behavior has been found in model 
calculations based on the exchange of spin fluctuations including processes not considered 
here. 47 ' 48 It is also possible that critical effects in a range of up to several degrees near the 
transition may lead to a divergence more rapid than in the usual mean field case. 4 ^ 

VI. Conclusions 

In this paper we have calculated A(0,T) and <r(0,T) within the framework of a BCS 
model in which the gap has d x 2_ y 2 symmetry, and both strong elastic impurity scattering 
and spin-fluctuation inelastic scattering processes are taken into account. We have sought 
to address a set of basic questions raised in the introduction. Here we summarize what we 
have learned. 

1) The microwave conductivity of the layered cuprates can be written in a Drude-like 



Here n qp (T) is the normal quasiparticle fluid density and the brackets denote the fre- 
quency average defined in Eq. (28). The inverse quasiparticle lifetime t _1 (u;, T) is the 
sum of the elastic impurity scattering rate and the inelastic spin-fluctuation scattering. 



form 




(34) 
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The form of Eq. (28) describes the transport properties of nodal quasiparticles which 
have a relaxation time t(uj, T) and a density of states N(uj). 

2) In the hydrodynamic limit O(r) <C 1, cr(T) = n qp (T)e 2 (T)/m. This is just the form that 
Bonn et al. used to extract a quasiparticle lifetime from their conductivity data. Here 
we have shown that (r) corresponds to an average over a frequency- and temperature- 
dependent lifetime. Figure 8 shows a plot of (r) _1 versus T for typical parameters. 

3) We find that for a d x 2_ y 2 gap, a(T — > 0) goes to a constant a ao = ne^/mnAQ indepen- 
dent of the impurity concentration (for small concentrations). 33 If we take r _1 (T c ) ~ T c 
from DC resistivity measurements, and 2Aq//cT c = 6, then a 00 /a(T c ) — 1/3% so that 
the limiting value of cr 00 is a about an order of magnitude smaller than cr(T c ). As 
the temperature increases, a(T) grows as T 2 . For T > T*, this can be understood as 
arising from the fact that both n qp (T) and (r) in the resonant scattering limit vary 
linearly with T. Note that we also find that for T < T*, a(T) - cr 00 varies as T 2 . If, 
in the pure limit T > T*, (r) were a constant, then cr(T) would increase linearly with 
T. However, this is not the case for the model we have considered. Both the fact that 
a(T — > 0) is independent of the impurity concentration and that <r(T) increases as T 2 
appear to be in disagreement with the presently available data. There is some evidence 
that the residual conductivity may be substantially lowered by reducing the density 
of twins in the crystal, 5 but the linear-T behavior remains a puzzle. Whether other 
scattering mechanisms can give rise to this behavior is not at present understood. The 
effect of particle-hole asymmetry is of particular interest in the context of our observa- 
tion that a constant relaxation time at low temperatures in pure samples is needed to 
produce a linear temperature dependence. The analytic properties of the self-energy 
of a particle- hole symmetric superconductor formally preclude such a result, however. 
An investigation of particle-hole asymmetry effects is in progress. 

4) At higher temperatures, inelastic scattering processes become important and give rise 
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to a scattering rate which increases initially as (T/T c ) 3 . As shown in Fig. 8, this leads 
to a minimum in (r) _1 at a particular value of T/T c . 

5) At higher microwave frequencies where ft(r) ~ 1, there is a crossover from the hydro- 
dynamic to the collisionless regime, and the relationship of cr(T, ft) to the quasiparticle 
lifetime involves an average of r(w,T)/(l + fiV(w,T)). In this regime, the conduc- 
tivity can exhibit a quasi-linear variation with T. We have shown in Fig. 10 how the 
temperature T p of the peak conductivity varies with ft along with cr(ft, T p )/a(T c ). We 
have also found that at higher microwave frequencies, quasiparticle screening leads 
to a reduction in X(T,ft). At a fixed temperature X(T,ft) can approach A(0, 0) as 
ft increases. We have used the full frequency dependence of A(T, O) and cr\{T, ft) in 
calculating the surface resistance R S {T, ft) shown in Fig. 16. 

6) In Section V, we explored the extent to which the d x 2_ y 2- wave plus scattering model 
can describe the surface impedance observed in YBa2Cu306 95 and its Zn-doped vari- 
ants. It appears (Figs. 11 and 12) that the temperature- and impurity- dependence of 
the penetration depth can be fit within the framework of this model. It will be interest- 
ing to compare the results for the O-dependence of A(T, ft) with experimental results 
which will soon be available. 50 The measured values of 01 (T, ft) shown in Fig. 14 for 
the pure and 0.15 % Zn samples appear to have a linear low-temperature variation in 
contrast to the T 2 variation predicted from the model. In addition, as noted, the lim- 
iting residual value of the conductivity obtained from the theory is smaller than that 
observed in many samples and is independent of the concentration of impurities. Nev- 
ertheless, as shown in Figs. 11 and 15, a simple <i-wave model plus scattering provides 
a reasonable overall fit to both the real and imaginary parts of the conductivity. One 
can ask whether alternative models such as an anisotropic s-wave pairing could provide 
similar fits to the data. In the absence of impurity scattering, the penetration depth 
and the low- frequency microwave conductivity cr(T) will both vary exponentially at 
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temperatures below the minimum gap value. In addition, if the minimum gap value is 
finite, cr(T — > 0) will vanish as exp — (A m j n /T). An extreme example of an anisotropic 
s-wave gap is given by taking for A the magnitude of the d x 2_ y 2 gap, Aq(T) \ cos20|. 
In this case, the results in the pure limit for A(T) are identical to the d x 2_ y 2 results. 
However, the addition of impurities can lead to a qualitatively different behavior for the 
anisotropic s-wave case. 51 As discussed in Section II, both £u n and A^, are renormalized 
by impurities in the s-wave case. In particular, potential scattering acts to average the 
gap over the Fermi surface, thus reducing the peak value of the gap and increasing the 
minimum value. Thus, even if one took the extreme anisotropic s-wave case in which 
the gap has nodes but does not change sign, impurities would lead to a finite effective 
gap and an exponential rather than T 2 crossover of the low-temperature dependence of 
both A(T) and cr(T). If "inert" defects like Zn impurities are found to have a magnetic 
character, 52 however, distinguishing s— and d— wave states becomes more difficult. 51 
Further measurements of the low-temperature dependence of the surface impedance 
in pure and impurity doped cuprates along with detailed comparisons with theoretical 
models are necessary to determine the symmetry of the pairing state. 
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Figure Captions 

1. Normalized low-T conductivity, a/aoo vs - the reduced temperature T/T c for microwave 
frequency O = 0. The solid lines correspond to resonant scattering, c = 0, T/T c = 
0.01, 0.003, 0.001, and dashed line corresponds to c = 0.3, F/T c = 0.01. 

2. Normalized low-T conductivity, (t/gqq vs. the reduced temperature T/T c in the Born 
limit, T N /T C = 0.01, 0/T c = 0, 0.001, 0.01. 

3. Normalized low-T conductivity, cr/aoo vs. the reduced temperature T/T c in the reso- 
nant limit, for c = 0, F/T c = 0.001, and 0/T c = 0, 0.0032, 0.01. 

4. Normalized conductivity, (j/ctqq vs. the reduced frequency 0/T c for T = 0, and T/T c = 
0.001,0.01,0.1. 

5. Normalized zero-temperature London penetration depth, A(T = 0)/Ao vs. the reduced 
scattering rate, T/T c q in the resonant scattering limit, c=0. 

6. Normalized London penetration depth, A(T)/Ao vs. the reduced temperature, T/T c 
for resonant scattering, F/T c = 0.0008, c = 0, and 0/T c = 0, 0.002, 0.018. 

7. Impurity relaxation rate 1/T c r(u) vs. the reduced frequency uj/Aq for Y/T c = 
0.01, 0.001 and c = (solid lines) and T/T c = 0.01, c = .2 (dashed line). 

8. Relaxation rate including inelastic scattering 1/T c (r) vs. the reduced temperature 
T/T c for T/T c = 0.0008, 0.009, 0.018, c = 0, A /T c = 3 (solid lines) and T/T c = 
0.0008, c = 0, Aq/T c = 4 (dashed line). 

9. Normalized conductivity including inelastic scattering, cr/aoo vs. the reduced tem- 
perature T/T c in the resonant limit, c = for 0/T c = 0.018, c = 0, and T/T c = 
0.0008,0.009,0.018. 

10. Reduced conductivity peak temperature, T p /T c vs. 0/T c for T/T c = 0.001, 0.01, c = 0, 
and Aq/T c = 3 (left axis); normalized peak conductivity a(T p , Q)/a(T c , 0) vs. 0/T c 
(right axis). 
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11. Comparison of d-wave penetration depth with penetration depth data on YBCO single 
crystals. 5 Normalized penetration depth, A(T)/Ao vs. the reduced temperature T/T c 
for r/T c = 0.018, 0.009, 0.0008 and c = 0. Data for pure YBCO crystal (circles), 0.15% 
Zn (diamonds), and 0.31% Zn (squares). 

12. Normalized T=0 normal fluid density 1 — (Aq/Aq) vs. the reduced coefficient of T 2 
term, Aq/(c2Aq) in the d x 2_ y 2 plus resonant scattering model. Each cluster of points 
represents one YBCO film from Ref. 45. 

13. Normalized theoretical conductivity a/ai(T c ) vs. the reduced temperature T/T c for 
impurity parameters T/T c = 0.0008 and c = 0, including inelastic scattering for 0/T c = 
0.002 and 0.018 (solid lines). Data points are normalized conductivities of YBCO 
single crystals from Ref. 5 for microwave frequencies 3.88 GHz (circles) and 34.8 GHz 
(triangles) . 

14. Normalized theoretical conductivity a/ai(T c ) vs. the reduced temperature T/T c for 
impurity parameters T/T c = 0.0008,0.009 and 0.018 with c = 0, including inelastic 
scattering for 0/T c = 0.018 (solid lines). Data points are normalized conductivities of 
YBCO single crystals from Ref. 5 for frequency 34.8 GHz, for samples nominally pure 
(circles), 0.15% Zn (triangles), and 0.31% Zn (squares). 

15. Effect of detwinning. Normalized theoretical conductivity a / a\(T c ) vs. the reduced 
temperature T/T c for impurity parameters T/T c = 0.0008 and c = 0, including inelastic 
scattering for 0/T c = 0.002 (solid line). Data points are normalized conductivities of 
detwinned YBCO single crystal from Ref. 5 for frequency 4.1 GHz. 

16. Normalized surface resistance, R s f R S (T C ) vs. the reduced temperature T/T c . Theory 
for 0/T c = 0.002 and impurity parameters T/T c = 0.0008, c = 0, including inelastic 
scattering, for Aq/T c = 3 (solid line) and Aq/T c = 4 (dashed line). Data from Ref. 5, 
3.88GHz, nominally pure YBCO crystal. 
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